J/ip dissociation by pions in QCD 
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We compute the J/tp n ^ charmed mesons cross section using QCD sum rules. This cross section 
is important to distinguish the suppression in the production of J/tp through the dissociation by 
comoving pions and through the formation of quark-gluon plasma. Our sum rules for the J/tp n — + 
D D* , D D* , D* D* and D D hadronic matrix elements are constructed by using vacuum-pion 
correlation functions, and we work up to twist-4. After doing a thermal average we get (a'^'^^^v) ~ 
0.2 - 0.4 mb at T = 150 MeV. 



In relativistic heavy ion collisions J/ip suppression has 
been recognized as an important tool to identify the 
possible phase transition to quark-gluon plasma (QGP). 
Matsui and Satz I*! predicted that in presence of quark- 
gluon plasma, binding of a cc pair into a J/ijj meson will 
be hindered, leading to the so called J/iJj suppression in 
heavy ion collisions. Over the years several experiments 
measured the J/tp yield in heavy ion collisions (for a re- 
view of data and interpretations see Ref.|3,l3l)- In brief, 
experimental data do show suppression 'i^. However, this 
could be attributed to more conventional J/ip absorption 
by comovers, not present in pA collisions [aSH- 
der to confirm that the suppression of J/ip comes from 
the presence of the QGP, it is necessary to understand 
better the J/ip dissociation mechanism by collision with 
comoving hadrons. 

Since there is no empirical information on J/ip ab- 
sorption cross sections by hadrons, theoretical models 
are needed to estimate their values. In general, differ- 
ent models apply to different energy regimes and one 
of the first estimates of the charmonium-hadron cross 
section uses short distance QCD |^ |^ 0|. However, 
the method is inapplicable at low energies, which is the 
regime of greatest interest for J/ tp collision with comov- 
ing hadrons. Besides, even in the high energy regime, 
nonperturbative effects may be important 11] and can in- 
crease significantly the value of the cross section. At the 
low energy regime one can use quark-interchange models 
[l^ or meson exchange models 14] . The results of 
the calculations for the charmonium-pion cross section 
based on these two approaches can differ by two orders 
of magnitude in the relevant energy range. Moreover, the 
rapid increase of the cross section near the threshold, ob- 
tained with these two models, is probably overestimated 
since these models do not respect chiral symmetry, as 
showed in [14 . There is also a calculation of the J/tp — n 
cross section lal based on the QCD sum rules (QCDSR) 
technique |l6l llTI ] , which is also valid at the low energy 
regime. The result for the cross section in ref. |l5j is in 
between the results in the quark-interchange models and 
meson exchange models. 

In this work we improve the calculation done in ref. [isj 
by considering sum rules based on a three-point function 



with a pion. We work up to twist-4, which allows us to 
study the convergence of the OPE expansion. Since the 
method of the QCDSR uses QCD explicitly, we believe 
that our work will improve the understanding of this im- 
portant topic. 

In the QCDSR approach, the short range perturba- 
tive QCD is extended by an OPE expansion of the cor- 
relator, giving a series in inverse powers of the squared 
momentum with Wilson coefficients. The convergence 
at low momentum is improved by using a Borel trans- 
form. The coefficients involve universal quark and gluon 
condensates. The quark-based calculation of a given cor- 
relator is equated to the same correlator, calculated us- 
ing hadronic degrees of freedom via a dispersion relation, 
giving sum rules from which a hadronic quantity can be 
estimated. 

Let us start with the vacuum-pion correlation function 
for the process J/tp tt — > Z) D: 
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with the currents given by — C7^c, jo = uij^c and 
jjj — cij^d. pi, P2, P3 and p4 are the four- momenta of the 
mesons tt, J/ip, D and D respectively. The advantage of 
this approach as compared with the 4-point calculation 
in ref. ,15j . is that we can consider more terms in the 
OPE expansion of the correlation function in Eq. and, 
therefore, we get a much richer sum rule. 

Following ref. [3], we can rewrite Eq. as: 



n„ 



d'^k 



Tr[S{p3 - k)j^ 



(27r)4 

X S{p3 - P2 - k)j5Daa{k,Plh5] 



(2) 



where S{p) is the free c-quark propagator, and Dab{k,p) 
denotes the quark-antiquark component with a pion, 
which can be separated into three pieces depending on 
the Dirac matrices involved (isf : 



Dab{k,p) = dab [ilbA + 7q75-B" + ^^aapO 



a/3] 
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The three invariant functions of {k,p): A, B" and 
C"^, are defined by the Fourier transform of the 
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vacuum-pion matrix elements: (0|d(x)i75M(0)|7r(pi)), 
(0|d»7"75u(0)k(pi)) and (0| J(x)a"''75^i(0)|7r(pi)) re- 
spectively. 

Using PCAC and working at the order O(ppPiy) we get 
up to twist-4 [llligt: 
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where tuq and 5^ are defined by (qD^q) = ra'^{qq)/2, 
mgsQ''^lpu\TT{p)) = iS^Up'' , with Ga^p = e„/5,,g-72 
and Gap = t^Gap- 

The additional contributions to the OPE comes from 
the diagrams where one gluon, emitted from the c-quark 
propagator, is combined with the quark- antiquark com- 
ponent. Taking the gluon stress tensor into the quark- 
antiquark component, one can write down the correlation 
function into the form 
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where we have defined 
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The c-quark propagator with one gluon attached is given 
by M gsGapS'^'^ik), and 

D^^ik,?) = j5<JpxEP^'"'{k,p)+l^^''^''Frgs{k,p) . 

the two functions 



Up to twist-4 and at order 0_ 
appearing above are given by 
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states contribution to the matrix element in Eq. Q . The 
hadronic amplitude is defined by the matrix element: 

iM {tp{p2,n)\ D{-p3) D{-pi) 7r(pi)) 

= i M^{pi,p2,P3,Pi) £2 • (8) 

The phenomenological side of the sum rule can be writ- 
ten as (for the part of the hadronic amplitude that will 
contribute to the cross section) [l5| : 
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where h. r. means higher resonances. The hadronic am- 
plitude can be parametrized as: 
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where A is the parameter that we will evaluate from the 
sum rules. 

Inserting the results in Eqs. 10} and Q into Eqs. Q 
and (jsj we can write a sum rule for the invariant structure 
appearing in Eq. Hl()(l . To improve the matching between 
the phenomenological and theoretical sides we follow the 
usual procedure and make a single Borel transformation 
to all the external momenta taken to be equal: —p2 = 
M^. We get, in the approximation 
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where we have transferred to the theoretical side the cou- 
plings of the currents with the mesons. The problem of 
doing a single Borel transformation is the fact that terms 
associated with the pole-continuum transitions are not 
suppressed ji^. In the present case we have two kinds 
of these transitions: double pole-continuum and single 
pole-continuum. In the limit of similar meson masses it 
is easy to show that the Borel behavior of the three-pole, 
double pole-continuum and single pole-continuum contri- 
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where /a^ is defined by the vacuum-pion matrix element 

mgs'J<.pi^G"^u\TT{p)) [i3. 

The phenomenological side of the correlation function, 
n^, is obtained by the consideration of J /tp, tt, D and D 



respectively. Therefore, we can single out the three- 
pole contribution from the others by introducing two 
parameters, A and -B, in the phenomenological side of 
the sum rule, which will account for the double pole- 
continuum and single pole-continuum contributions re- 
spectively |ilEill2a|. 

The parameter values used in all calculations are 
TOc = 1.37 GeV, = 140 MeV, mo = 1.87 GeV, 
mo' = 2.01 GeV, m^ = 3.097 GeV, = 131.5 MeV, 
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FIG. 1: Amplitude of the process J/t/j tt D D as a, function 
of the Borel mass. The circles, squares and diamonds give the 
twist-2, 3 and 4 contributions to the sum rule. The triangles 
give the result from Eq. IIH . The solid line give the fit to the 
QCDSR results. 



240 MeV, 
0.2 GeV^ 



= 270 MeV, fo = 170 MeV, fo- = 
(qq) = -(0.23)3 GeV^ TO§ = 0.8 GeV^, 
= 0.0035 GeV^ 0- 

In Fig. 1 we show the QCD sum rule results for 
A + AM'^ + BM'^ as a function of M^. The circles, 
squares and diamonds give the twist-2, 3 and 4 contribu- 
tions respectively. The triangles give the final QCDSR 
results. We see that the twist-3 and 4 contributions are 
small as compared with the twist-2 contribution show- 
ing a "convergence" of the OPE expansion. The tri- 
angles follow almost a straight line in the Borel region 
6 < < 16 GeV^, indicating that the single pole- 
continuum transitions contribution is small. The value 
of the amplitude A is obtained by the extrapolation of 
the fit to = |2(T[ElIEi|. Fitting the QCD sum rule 
results to a quadratic form we get A ~ 11.4 GeV"^ As 
expected, in our approach A is just a number and all de- 
pendence of A^^ on particle momenta is contained in the 
Dirac structure. This is a consequence of our low energy 
approximation. 

Instead of using the experimental values for the me- 
son decay constants, it is also possible to use the re- 
spective sum rules, as done in [l5j |. The behavior of 
the results does not change significantly, leading only to 
a change in the value of the amplitude. Using the re- 
spective sum rules for the meson decay constants we get 
A ~ 14.9 GeV"^ We wiU use these two procedures to 
estimate the errors in our calculation. Our results agrees 
completely with the value obtained in |l5l |. 

The calculation of the sum rules for the processes 
J/ip TT ^ D D* and J/tp tt ^ D* D* can be done in 
a similar way. One has only to change the currents in 
Eq. by the appropriate ones. The hadronic ampli- 



tudes for these two processes can be written in terms of 
many different structures. In terms of the structures that 
will contribute to the cross section we can write 
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for the process J/ip tt — > D* D*: 
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£p.pai3P2i'PiP2 and 



£tiua0Pip)P2P3 + <^pvai3P2pPipl- In principle aU the 
independent structures appearing in H^^p and J^^p 
would have independent parameters A^. However, since 
in our approach we get exactly the same sum rules for 
all of them, we decided to group them with the same 
parameters. 

The exp ressions for all 20 sum rules will be given else- 
where 1231 . At this point it is important to stress that the 
sum rule for the process J / tp ir ^ D D is not particular, 
in general all the other sum rules are similar and contain 
twist-2, twist-3 and twist-4 contributions corresponding 
to the first, second, and third terms inside the brackets 
in the right hand side of Eq. Hll|l . Only the sum rules for 
Af^ ^ up to A{4 ^ get only twist-4 contributions, and 
give results compatible with zero. It is also interesting to 
notice that if we consider only the twist-2 contributions 
we recover the sum rules obtained in ref. |l5l |. 

The results for all other sum rules show a similar be- 
havior and the amplitude can be extracted by the ex- 
trapolation of the fit to = 0. The values for all the 
parameters are given in In Eq. (|12|l the structures 
multiplying A4 and A5 break chiral symmetry 1^ and, 
therefore, will be neglected. 

Having the QCD sum rule results for the amplitude of 
the three processes J/ip n ^ D D*, D D, D* D* , we can 
evaluate the cross section. In Fig. 2 we show the cross 
section for the J/t/j — tt dissociation. The shaded area 
give an evaluation of the uncertainties in our calculation 
obtained with the two procedures described above. It is 
important to keep in mind that, since our sum rule was 
derived in the limit pi « p2,P3,PA, we can not extend 
our results to large values of ^/s. 
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FIG. 2: Total J/t/j tt dissociation cross sections of the pro- 
cesses J/ip -K ^ D D* + D D* + D D + D* D*. The shaded 
area give an evaluation of the uncertainties in our calculation. 

In a hadron gas, pions collide with the J/ip at different 
energies. The momentum distribution of thermal pions in 
a hadron gas depends on the effective temperature T with 
an approximate Bose-Einstein distribution. Therefore, 
the relevant quantity is not the value of the cross section 
at a given energy, but which is the product of 

the dissociation cross section and the relative velocity 
averaged over the energies of the pions. 
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FIG. 3: Thermal average of J/^ dissociation cross section by 
pions as a function of temperature T. The shaded area give 
an evaluation of the uncertainties in our calculation. 

As shown in Fig. 3, (a'^'^^^v) increases with the tem- 
perature. Since the J/tp dissociation by a pion requires 
energetic pions to overcome the energy threshold, it has a 



small thermal average at low temperatures. The shaded 
area in Fig. 3 give an evaluation of the uncertainties in 
our calculation due to the two procedures used to extract 
the hadronic amplitudes. 

In conclusion, we have studied the J/^ dissociation 
cross section by pions using the QCDSR technique, based 
on a three-point function using vacuum-pion correlation 
functions. We have estimated the hadronic amplitudes 
by working up to twist-4 in the limit pi << P2,P3,P4- 
Our results are in agreement with the former QCDSR 
calculation, done with a four-point function at the pion 
pole Our results for the cross section as a function 
of ^/s are smaller than the results using meson-exchange 
models (without form factors), but larger than the cal- 
culation based on quark-exchange models. 

The dominant contribution to the hadronic amplitudes 
comes from the twist-2 operator, or equivalently, from the 
quark condensate. As we know that the quark conden- 
sate is stronger in the vacuum and weaker in the interior 
of hadrons, we can conclude that the charmonium "sees" 
and interacts with the surface of the pions, where there 
is a "halo" of condensates. This is way the cross sec- 
tion can be larger than the geometric value. In our ap- 
proach the continuous growth of the cross section comes 
from the growth of the phase space, as in the effective 
Lagrangian calculations. In the short distance QCD cal- 
culation 0, ^3 the cross section also grows with y/s, but 
the growth there is considerably smaller because the rise 
in the gluon density cannot completely compensate the 
fall of the partonic cross section Jj| . 

The thermal average of the J/ip — tt dissociation cross 
section increases with the temperature and at T = 
150 MeV we get {a'^^l'^v) ~ 0.2-0.4 mb which is smaller 
than the values used in phenomenological studies of J ji]} 
absorption by comoving hadrons in relativistic heavy ion 
collisions. 

The same approach used here could be ap plie d to cal- 
culate the T TT cross section. In a recent work [iJJ the cttit 
was computed using short distance QCD. Since we expect 
the non-perturbative corrections to be less important for 
heavier systems, the differences between short distance 
QCD and QCD sum rules should be smaller for the Ttt 
system, and a systematic comparison between the two 
approaches becomes possible. We will address this point 
in the future. 

Another possible extension of this work is the calcu- 
lation of the Xc TT and r]c n cross sections. This can be 
done by replacing the current in Eq. by the corre- 
sponding Xc and rjc currents. Unfortunately since ip' and 
J/ip have the same quantum numbers and are, therefore, 
described by the same current, it is not possible, in this 
approach, to estimate the ip' tt cross section. The ■0' con- 
tribution to the present sum rule calculation is inside the 
parameters A and B, in Eq. and cannot be separated 
from the other higher mass states contributions. 
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